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ABSTRACT 

SBS 0335-052 is a blue compact dwarf galaxy (BCD) with one of the lowest known metallicities, 
Z~Z0/41, making it a local example of how primordial starburst galaxies and their precursors might 
appear. A spectrum obtained with the Infrared Spectrograph (IRS 1 ) on the Spitzer Space Telescope 
clearly shows silicate absorption features, emission lines of [SIV] and [Nelll], and puts strong upper 
limits on the PAH emission features. The observed low resolution spectrum (R~90) extends from 5.3 
to 35 /xm and peaks at ^28 /im. The spectrum is compared to IRS observations of the prototypical 
starburst nucleus NGC7714. SBS 0335-052 is quite unlike normal starburst galaxies, which show 
strong PAH bands, low ionization emission lines, and a continuum peak near 80 /im. The continuum 
difference for A > 30 /im implies a substantial reduction in the mass of cold dust. If the spectrum of 
this very low metallicity galaxy is representative of star forming galaxies at higher redshifts, it may 
be difficult to distinguish them from AGNs which also show relatively featureless flat spectra in the 
mid-IR. 
Subject headings: dust, extinction — galaxies: individual (SBS 0335-052) — galaxies: starburst 



1. INTRODUCTION 

By virtue of their low metallicities, small size and 
low optical luminosity (Mb >-18), Blue Compact Dwarf 
galaxies (BCDs) may be local examples o f the building- 
blocks for the earliest galaxies l|Reesll9 98') . A large num- 
ber of BCDs have been identified based on their blue 
colors and very low metallicities as determined by opti- 
cal spectroscopy, which implies that star formation has 
begun only recently (see IKunth fc Ostlinl 120001 for a 
review). Their blue color arises from one or more in- 
tense regions of active star formation that often appear 
nearly devoid of d ust. The first identified member of this 
class was I Zw 18 l|Searle fe Sargenall9 72'). which still re- 
mains the most metal poor member of the class with 
Z=Z0/5O. Another object wit h similarly low me tallic- 
ity Z=Z©/41, is SBS 0335-052 ijlzotov et al.l H997) at a 
distance 5 of 57.6 Mpc. In the case of SBS 0335-052 six 
regions of massive star formation, five visible and one ob- 
scured, have been identified, with ages less than 25Myr 
and a total luminosity of ~ 10 9 L^. All six are withi n 
a region of ~2" or — 500pc in size l|Thuan et alJll997j) . 
but roughly 75% of the tot al luminosity comes as m id- 
infrared (mid-IR) radiation (Plante & Sauvage 2002)! A 
few similar cases of obscured super star clusters which 
contribute a considerable fraction of the total luminos- 
ity of the host galaxy in the mid-IR have been observed 
so far. These include the Antennae interacting galaxies 
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(Mir abel et all 11 9981) and the Wolf-Rayet galaxy He 2- 
10 (|Vacca et al.ll2002|) . The frequency and implications 
of this phenomenon at higher redshifts (z>0.1) will re- 
quire higher sensitivity observations of the mid-IR emis- 
sion than have been available to date. SBS 0335-052 
was known to have an infrared-emitting dust continuum 
(Thuan et al. 1999), and its dust properties have already 
been studied quite extensively with both the Infrared 
Space Observatory (ISO) and from ground based tele- 
scopes [|Dale et all fcOOl: Pl ante fc Sauvagdl2002|) . How- 
ever, there is still disagreement regarding whet her the 
emitting dust is optically thick ( T huan et alJ 11999(1 or 
optically thin (Dal e et al1l2001|) . This impacts the ques- 
tion of how enshrouded the first generations of stars can 
be and whether or not there may exist a population of 
"optically-quiet" objects which are visible in the infrared 
but not in the optical. 

We chose SBS 0335-052, along with approximately 20 
other well studied BCDs, to investigate the characteris- 
tics of star formation a t very low metallicities using th e 
IRS l|Houck et al.ll2004]) on Spitzer IjWerner et alJ l20041. 
Similar low metallicity objects may be detected at much 
higher luminosity and much greater redshifts in Spitzer 
discovery surveys. 

2. OBSERVATIONS 

SBS 0335-052 was observed using both IRS low resolu- 
tion modules. The spectrum extends from 5.3 to 35 /im 
and was obtained on 6 February 2004. The red peak-up 
camera in medium-accuracy mode was used to locate the 
mid-IR centroid of the source and move it to the center 
of the spectrograph slits. The total integration time was 
42 minutes with 28 minutes for the Short-Low module 
and 14 minutes for the Long-Low module. The total 
elapsed time including the telescope slew, target acqui- 
sition, settling, array conditioning, and integration was 
61.3 minutes. 

The basic processing of the data, such as ramp fit- 
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ting, dark sky subtraction, removal of cosmic rays, droop 
and linearity correction, wavelength calibration, etc, was 
performed using the IRS pipeline at the Spitzer Science 
Center (version S9.1). The resulting spectral images were 
sky-subtracted and a one-dimensional spectrum then ex- 
tracted. The peak-up images were also used to derive 
a photometric point at 22 /im (filter bandwidth 18.5- 
26.0 /mi). As described in detail in chapter 7 of the 
Spitzer Observer's Manual 6 , during IRS peak up we ob- 
tain 6 images of the science target. The on source time 
for the 22 /im peak- up images of our target was 6x8=48 
seconds and each image was created by reading the ar- 
ray in double-correlated sample mode. We processed the 
data on the ground to remove cosmic rays and the resid- 
ual noise of the electronics. The resulting image had 
a prominent diffraction ring and was indistinguishable 
from the image of a point source. The conversion to flux 
density was based on a number of calibration stars for 
which peak-up imag es, IRS spectra, a nd reliable tem- 
plates are available IjCohen et alJl2003|) . We find that 
the 22 /im flux density of SBS 0335-052 is 70±llmJy. 

3. RESULTS 

3.1. Mid-IR Spectral Properties 

Figure 1 shows the 5.3-35 /mi spectrum of SBS 0335- 
052 as observed by the IRS. Our data are in good agree- 
ment with the ov erall shape and inten sity of the 5-15 /mi 
ISO spectrum of iThuan et alJ y,999), but our signal to 
noise is at least a factor of 10 higher. This enables us 
for the first time to directly detect a few mid-IR ionic 
lines, while placing strong upper limits on others. The 
9.7 /zm silicate absorption feature is clearly evident in our 
spectrum, and the 18 /im feature is probably present. 
Using the silicate a bsorption profile meas ured towards 
the Galactic center (jChiar fc Tielensll2004D we find that 
for a screen model Ag.7 M , n = 0.49 mag, assuming a black- 
body background source. More importantly, our spec- 
trum indicates that the emission from the galaxy peaks 
at ^28 /mi which, as we discuss in the following sec- 
tion, has important consequences in estimating the dust 
mass and grain size distribution. In Figure 1 we also 
present a scaled version of the IRS low-resolution spec- 
trum for the protot ype starburst nucleus of NGC 7714 
IjBrandl et al J 120041) . normalizing its flux to the corre- 
sponding flux of SBS 0335-052 at 14 /im (the actual flux 
for NGC 7714 is ~9.5 t i mes la rger than shown). As also 
noted by IThuan et alJ l|1999|) , a striking difference be- 
tween the spectrum of SBS 0335-052 and that of a more 
typical starburst is the absence of strong polycyclic aro- 
matic hydrocarbon (PAH) features and low excitation 
ionic lines. In starburst galaxies, emission from PAHs 
is thought to originate from the photodissociation en- 
velopes bordering the H II regions produced by the ion- 
izing starburst. How these features can be absent in a 
low-metallicity starburst is an important astrophysical 
question. One possibility is that the absence of PAHs 
is due to low abundance of carbon and/or nucleating 
grains; another possibility is that PAHs are quickly de- 
stroyed. Table 1 shows that the observed line ratios of 
log([SiV]/[SIII])>0.48, and log([NeIII]/[NeII])>0.69 are 
similar to the most extreme example of ultra-compact 
H II regions in our Galaxy. The ratios indicate that the 

6 http://ssc.spitzer.caltech.edu/documcnts/som/ 
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Fig. 1.— The IRS low resolution spectrum of SBS 0335-052 is 
presented (solid line) alo n g with the spectrum of NGC 7714 (dotted 
line) from IBrandl et alJ 120041) . The spectrum of NGC 7714 has 
been divided by 9.53 so that its 14 fim flux density matches that of 
SBS 0335-052. Note the complete lack of PAHs in the spectrum of 
SBS 0335-052 as well as the spectrum peak at ~30 fira (see Section 
3.1). The solid circle is our 22 fira peak-up photometric point, 
while the diamond corre sponds to the 21 fim Gemini point from 
IPlante fc Sauvage! 120021) . 



radiation field is extremely hard and corresponds to an 
effective stellar temperat ure of T e f f > 4 x 10 4 K, a ssum- 
ing solar abundance (see iMartin- Hernan dez et alJ l2002). 
This would suggest that the absence of PAHs results from 
their destruction by the hard UV photon s and strong, 
wind s produced by the massive stars (i.e. lAllain et alJ 
1996). Such a scenario is quite likely since in low metal- 
licity systems the attenuation of UV photons is small 
and consequently their mean free path in the interstellar 
medium can be larger and photodissociation of PAH may 
occur over considerably larger scales than those seen in 
typical Galactic H II regions. Whatever the PAH life- 
cycle is in SBS 0335-052, the primary result is the well- 
defined absence of these features in the mid-IR. 

Our spectrum displays strong [SIV]A10.51/mi, and 
[NeIII]A15.55/mi lines, while [SIII], [Nell], and H 2 0-0 
S(3) may also be present (see Table 1). The observed 
(not corrected for extinction) [Nelll] and [SIV] fluxes 
in conjunction with the free-free radio e mission have 
been used to derive the ionic abundances l)Osterbrockl 
1989). This method is a priori more robust than us- 
ing hydrogen recombination lines because it is relatively 
free of complicating factors: extinction, uncertainties 
about electron temperatures and densities as well as un- 
derlying stellar absorption and wi nd associated emis - 
sion features. Using the results of iHunt et alJ (2004), 
which suggest a free-free emission of S^=0.27mJy at 
5 GHz and electron density of ~20 00cm~ 3 , as well as 
the helium abundance as given by [izotov et al. (1997) 
we find that Ne++/H+=1.05xl0- 5 which is ~9% of 
the solar value and S 3 +/H+=4.3xlQ- 7 or ~2% of so - 
lar IjAnders fc Grevessd Il989t iGrevesse fc Noels! Il993|) . 
These estimates are lower bounds on the total Ne and 
S abundances because only a single ionic state has been 
measured and extinction corrections have not been made. 
Had we corrected for extinction based on the observed 
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TABLE 1 
Mid-IR Line Fluxes of SBS 0335-052 



ion 


X(fim) 


Flux 


(xlO- 17 W m" 2 ) 


[SIV] 


10.51 




1.62±0.09 


[SIII] 


18.71 




<0.54 a 


[Nell] 


12.81 




<0.28 a 


[Nelll] 


15.55 




1.40±0.08 


H 2 0-0 S(3) 


9.67 




<0.37 a 



SBS0335-052 (IRS) 

Tokeuchi et al. 2003 (model) 



a Thc upper limits are 3-cr. 



9.7/im silicate absorption, the corresponding Aio.5i^ m = 
0.33 mag and Ai5.56 Alm = 0.072 mag would imply that the 
abundances could be higher by a factor of 1.36 and 
1.07 respectively These abundances are higher than the 
Zq/41 metallicity that is derived for oxygen from the op- 
tical observations and imply that the region probed by 
our spectrum has polluted its environment by mass loss 
from supernovae (SNe) and Wolf-Rayet stars. 

However, there are a number of issues associated with 
the above mentioned abundances which must be con- 
sidered. There is an inconsistency between the radio 
measurements and the fluxes of the hydrogen recom- 
bination lines quoted in th e literature. If we use the 
5 GHz free- free emission of IHunt et al. ( 2004) to pre- 
dict the corresponding extinction corrected H/3 flux we 
find a value of 6.1xl0 _14 ergcm~ 2 s _1 which is lower 
than the observed (not correcte d for extinct i on) v alue of 
S.SxlO-^ergcm^s- 1 bv llzotov k Thua n (1998). The 
predicted Bra fl ux is also a factor o f ~2-5 less than the 
one observed and lHunt et al.1 l)2004[) advocate that winds 
could account for some of the Bra flux. It remains un- 
clear though why these winds would not affect the other 
recombination lines. Even if we were to consider possi- 
ble errors in the estimate of the thermal radio emission 
the puzzle remains. Clearly a fraction of the observed IR 
line fluxes originates from the dust-free "optical" region 
of the galaxy. To estimate this we can use the optical 
[NeIII]A3868A measurements an d derived physical c on- 
ditions of the interstellar gas by llzotov et al.l l)1997fl to 
predict the co rresponding j NcIlT|A15.55/mi flux. For a 
T e =20,000K (|Izotov et al.l 119971 derives 19,200 K for 
O III) we find that the [NeIII]A15.55/zm flux predicted 
by the optical measurements is 0.75xl0 _14 ergcm~ 2 s" 1 , 
~2 times lower than what is observed. This would sug- 
gest that indeed there is a star forming region of the 
galaxy which is not accessible in the optical. The result 
though depends strongly on the electron temperature. 
If T e =15,000K then the predicted mid-IR flux is simi- 
lar to the one observed. However, this would place the 
Ne ++ region at the same temperature as the + region; 
an unlikely situation. The 9.7/im silicate feature clearly 
indicates that there is a powerful source of luminosity 
hidden behind a screen of at least A. 
or Ay ~15mag. 



9.7^111 ~ 0.5 mag 
Clearly all of these uncertainties can 
strongly affect any abundance calculations. 

Another major difference between the two spectra 
shown in Figure 1 is that in SBS 0335-052 the con- 
tinuum shortward of ~15 /im appears similar to that 
which underlies the emission features in NGC 7714. How- 
ever, these continua depart dramatically at longer wave- 
lengths. The continuum of NGC 7714 increases rapidly at 
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Fig. 2.— The IR S spectrum of SBS 0335-052 (solid line) 
along with the model of Takeuchi ct al. (2003) (dotted line). Note 
how the model deviates past 16 /im. Based on the observed 5— 
35 fim SED of the galaxy we can extrapolate to longer wavelengths 
(dashed line) using a an offset power-law peaking at 28 (im with an 
asymptotic shape, fviy — * 0) ~ u 1 ^. The difference between the 
observed SED and the model has only a small effect in the infrared 
luminosity, but it substantially reduces the amount of cold dust in 
the galaxy. We also show estimates derived from the radio data of 
the average flux density over the far-infrared range and the upper 
limit at 60 (im (box and triangle, respectively; see Section 3). 



longer wavelengths because of a relatively more massive 
cool dust component which characterizes many luminous 
infrared galaxies. Conversely, the spectrum of SBS 0335- 
052 peaks at 28 /im. In Figure 2 we have drawn an offset 



power-law, f v 



,1.3 



, to the IRS spectrum to extrapolate 



to longer wavelengths. Using the non thermal 1.46GHz 
flux density of the galaxy measured by Hunt ct al. ( 2004) 
and the radio to far-infrared correlation, we predict that 
the far-infrared luminosity is F/i r (43-123Ltm)=2xl0~ 15 
Wm~ 2 . The corresponding average flux density over this 
wavelength range is ^44mJy. Even if the entire far- 
infrared luminosity originates from the 60 /xm band, the 
cor responding flux density w ould be 58mJy, well below 
the iPlante &: Sauvaed l)2002|) 112mJy measurement at 
65 /im. However, the 58mJy estimate is consistent with 
the IRS spectrum which clearly decreases at wavelengths 
longer than 30 /im. 

These two major differences, no PAH features and rel- 
atively flat continuum, provide an initial indication that 
the PAH criterion for identifying starbursts may not ap- 
ply to systems of low metallicity. Mid-IR spectra of Ac- 
tive Galactic Nuclei (AGN) often display ionic lines and 
also lack PAH emission features, but the continua are 
typically even flatter between 5-20 /im than SBS 0335- 
052 (see ISO observations of Seyfert Is by Clavel et al. 
2000, or the nuclear spectrum of NGC 1068 by Sturm et 
al. 2000 and Le Floc'h et al. 2001 as well as Peeters et al. 
2004). If metal poor galaxies such as SBS 0335-052 can 
be present at z~0, it is likely that similar systems may 
exist not only as high-z primordial galaxies, but also at 
intermediate redshifts 0.5 <z< 1.5 where moderate lumi- 
nosity galaxies (Ltr, < 10 11 L©) are accessible with the 
deep infrared surveys performed by Spitzer. Far-infrared 
and submillimetcr photometric methods for estimating 
redshifts based on spectral energy distributions (SEDs) 
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of nearby metal rich systems, with continua peaking at 
—100 /mi, would fail to identify metal poor systems with 
SEDs similar to that of SBS 0335-052, which peaks at 
-28 /mi. 

3.2. IR luminosity and Dust Mass 

The luminosity of a dust obscured source is the inte- 
grated infrared flux corrected for the object's distance. 
However, determining the total dust mass depends cru- 
cially on the dust' s temperature , distr ibution, and the 
optical properties. iThuan et alJ IJ1999J) concluded from 
their 5-15 /tm ISO spectrum that the starburst respon- 
sible for heating the dust which produces the mid-IR 
continuum in SBS 0335-052 is heavily obscured, with an 
optical extinction Ay— 20 mag. They also modeled the 
dust continuum as being represented by a single temper- 
ature of about 260K and an emissivity, e ~ v 15 . On 
this basis they suggested that the young star clusters 
visible optically were independent of the obscured clus- 



ters. These gener al conclusions have been su pport ed by 
Hunt et alJ MM. 12001 . iPlante fc Sauvagel IJ200l . and 



Takeuchi et al. (2003). If correct, this interpretation im- 



plies that primordial starbursts of very low metallicity 
can nevertheless surround themselves with sufficient dust 
to be completely hidde n optically. A ver y different con- 
clusion was reached bv lDale et alJ l)2001|) . using ground- 
based mid-IR imagery; they conclude that the optical ex- 
tinction to the starburst which heats the dust is less than 
one magnitude and that a two component dust model is 
needed, with temperatures of about 80K and 210K. Both 
components are optically thin, and the bolometric lumi- 
nosity is dominated by the cooler component. 

All previously published models extending beyond 
20 /mi had their assumed spectra anchored to the IRAS 
and ISO limits at 6 and 100 ^m, as well as to the ISO 
detec tion at 65 /mi l)Dale et alJl2001i IPlante fc Sauvagg 
2002). The shape and intensity of the IRS spectrum 
clearly demonstrates that those far-infrared values do not 
accurately reflect the true shape of the spectrum, which 
is inconsistent with a flux density above — 70mJy any- 
where in the 30-100 /tm range (see Figure 2). We have 
fit both the new IRS spectrum and the old ISO spec- 
trum with a simple two-temperature model with emis- 



sivity of e 



,1.5 



In both cases good fits are obtained 



with dust temperatures of —65 and — 150K. While the 
— 150K component is the same for both models, the 
required dust mass for the cool component is 4 times 
less for the Spitzer spectrum. Scaling to the predicted 
100 /im flux, we estimate the mass of the cool dust to be 
6xl0 3 M for the spectrum base d on the ISO data and 
1.5xl0 3 M Q for the IRS spectrum l)Hildebrandll983HCoxl 
2000). These masses are considerabl y below the masses 
reported previously (~ 10 5 M^, see IThuan et alJll999t 
iHunt et aill2001i IPlante k. Sauvagdl2002t) . Qualitatively 
this is because the cold component in the Spitzer model 



is warmer and produces less flux than the cold compo- 
nent in the ISO models. Dramatic changes in the dust 
properties could increase the dust mass, but the materi- 
als used in the cited models do not differ greatly from the 
stand ard. Several models, including the one developed 
by Takeuchi et al. (2003), derive a large fraction of the 
dust components from the SNe ejecta produced in the 
enshrouded super star clusters. However, the amount 
of dust produced by SNe in lo w metallicity systems i s 
predicted to be rather low (i.e. iTodini fc Ferraral l2001). 
Furthermore, in SNe remnants of our Ga laxy most of 
the d ust appears to be very cold (T~ 18K. lDunne et alJ 
2003) . The fact that the Spitzer data suggest that much 
less cool dust is required eases the problem of accounting 
for the in situ dust generation. 

4. CONCLUSIONS 

The very low metallicity Blue Compact Dwarf galaxy 
SBS 0335-052 is shown to have a very unusual spectrum 
which is quite different from the spectrum of typical star- 
burst galaxies — the flux density, f„, peaks at — 28 /tm 
while the luminosity, vi u , peaks at — 20 /im. There are 
no detectable PAH emission features. The spectrum is 
characterized by a warm (~ 150K) dust component with a 
more massive cool (-65K) dusty envelope. However, the 
mass of the cool region is far less than what had been pre- 
viously estimated. Silicate absorption features are clearly 
present (Ag^m > 0.49 mag). Preliminary analysis of the 
infrared ionic lines suggests that the young central clus- 
ter may have already polluted the cocoon enshrouding it, 
and its metallicity could be higher than that determined 
from optical observations alone. Upcoming mid-IR spec- 
troscopy using the high resolution modules of IRS will 
provide higher sensitivity and more accurate line fluxes. 

Theoretical models have indicated that the dust re- 
sponsible for the mid-IR emission can be explained by 
ejecta from type 2 SNe and Wolf-Rayet stars in the 
buried super star cluster. However, there are no com- 
plete models to date that fit the mid-IR spectrum and 
include dynamics of the dust generation and entrapment 
mechanisms. 
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